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ABSTRACT 



We present the first chemical analysis of stars on the double subgiant branch (SGB) of the globular cluster NGC 1851. We obtained 
48 Magellan IMACS spectra of subgiants and fainter stars covering the spectral region between 3650-6750A, to derive C and N 
abundances from the spectral features at 4300A (G-band) and at ~ 3883A (CN). We added to our sample ~ 45 unvevolved stars 
previously observed with FORS2 at the VLT. These two datasets were homogeneously reduced and analyzed. We derived abundances 
of C and N for a total of 64 stars and found considerable star-to-star variations in both [C/H] and [N/H] at all luminosities extending 
to the red giant branch (RGB) base (V~ 18.9). These abundances appear to be strongly anticorrelated, as would be expected from 
the CN-cycle emichment, but we did not detect any bimodality in the C or N content. We used HST and ground-based photometry to 
select two groups of faint- and bright-SGB stars from the visual and Stromgren color-magnitude diagrams. Significant variations in the 
carbon and nitrogen abundances are present among stars of each group, which indicates that each SGB hosts multiple subgenerations 
of stars. Bright- and faint-SGB stars differ in the total C-l-N content, where the fainter SGB have about 2.5 times the C-l-N content of 
the brighter ones. Coupling our results with literature photometric data and abundance determinations from high-resolution studies, 
we identify the fainter SGB with the red-RGB population, which also should be richer on average in Ba and other j-process elements, 
as well as in Na and N, when compared to brighter SGB and the blue-RGB population. 

Key words, stars: abundances - stars: sub giant branch -GCs: individual (NGC 1851)- C-M diagrams 



1. Introduction 

In the last years, both spectroscopic and photometric evidence 
has shown that globular clusters (GCs) cannot any longer be 
considered to be a simple stellar population, because self- 
enrichment is a common feature among them. While the detec- 
tion of significant scatter in iron and/or «-capture el ements asso- 
ciated to i-process is limited to a few clusters (e.g..' Hesser et al.l 
Ll982; Yong & Grundahl 2008; Marino et al. 2009), star-to-star 
variations in abundances of the light elements (C, N, O, Na, Mg, 
Al), have been observed in stars of al l evolutionary phases in the 
majority of GCs (e.g.,^artel l 2'oTltlPancino et al.boiOtlKrafi 
T979tlRamu-ez & Cohenl2003l;lKavser et al.l2008t Carretta et al 



2009) 



These variations manifest themselves through correlations 
and anticorrelations, the signature of high-temperature proton 



* Based on observations taken with the 6.5 meter Magellan Telescope 
at Las Campanas Observatory, Chile and with the Very Large Telescope 
at ESO La Silla Paranal Observatory, Chile, under programme ID 68.D- 
0510. 



fusion cycles that have processed C and O into N, Ne into Na 
and Mg to Al. The temperatures required to convert Ne into 
Na are on the order of T ~ xlO^ K and these are not reached 
in <0.8Mo GC dwarfs, which, furthermore do not posses a 
deep convective layer The generally accepted explanation is that 
these stars were born with the observed CNONa abundance pat- 
terns. Intermediate-mass asymptotic giant branch (AGB) stars, 
fast rotating massive stars, or massive interactive binaries have 
been proposed as sources of the necessary pollution of the intra 
cluster medium (I CM) before the second-gene r ation stars were 
formed (see, e.g. 'D'Anto na & Ventural 120071; iDecressin et al.l 
I2007t Ide Mink et al...2009,) . 

Star-to-star variations in light- and alpha-element abun- 
dances, age, and metallicity can determine multimodal or broad 
sequences in the CMD. Complex structures along the main se- 
quence (MS), subgiant branch (SGB), red giant branch (RGB) 
or horizon tal branch (HB) with i n some galactic or extra- galactic 



GCs ( e g-. iPancino et al 



20001; iBedin e t al. 2004; Solli ma et all 



120071: iPiotto et al.ll2007L Marino et al.ll2008[ iMilone et al.ll2010l: 
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iLardo et ani201 it to name a few), unambiguously indicate that 
GCs host two or more generations of stars. 

NGC 1851 is one of the most intriguing clusters with mul- 
tiple stellar populations. It exhibits a do uble SGB, with the 
faint component made of ~35% of stars dMilone et alJ l2008l 
hereafter M08). If age is the sole cause, then the SGB split 
is consistent with two stellar groups with an age difference of 
~1 Gyr. As an alternative, the two SGBs could be nearly co- 
eval but with a different C+N+O content dCassisi etal J 120081: 
I Ventura et aDl2009h . The HB is also bimodal, with about -35% 
of HB stars on the blue side of the instability strip. Both the 
HB and th e MS morph ology le ave no room f or strong helium 
variations dSalaris et al..,2008: Antona et aPHoogi). Nor is the 
RGB consistent with a simpl e stellar pop u lationlGrundahl et al 



A(RA) [arcmin] 




\T99^ ICalamida et alJ I2U07I) . iLeeetalJ (l2009h and iHanetal 



( 2009h pointed out two distinct RGB evolutionary sequences, us- 
ing Stromgren Ca uvby photometry, and proposed that the split 
might be attributed to differences in calcium abundance. 

Many spectroscopic studies have bee n dedicated to RGB 
stars in NGC 1851. Almost 30 years ago. iHesser et all (Il982h 
noticed three out of eight bright-RGB stars with anomalously 
strong C N bands and with enhanced Sr and Ba lines. More 
recently, lYong & Grundahll (l2008h analyzed UVES spectra of 
eight giants. Their analysis revealed that star-to-star abundance 
variations of O, Na, and Al with a clear anticorrelation of O and 
Na also exist in this cluster, and the amplitude of these variations 
is comparable with those found in clusters of the same metal- 
licity. More interestingly, they argued for the presence of two 
different populations in this cluster, characterized by significant 
differences in the the li ght ^-proce ss element Zr and the heavy 
i-process element La. lYong et aP (l2009l) and Yong et al. (in 
preparation) found a wide spread in the abundan ce sum C-hN-hO 
(whil e a constant sum of Ch-N+O was derived bv lVillanova et al.l 
I2OIOI from the abundance analysis of 15 RGB stars); with the 
CNO-rich stars being also enhanced in Zr and La. Yong and col- 
laborators associated the group of CNO-rich i-rich stars to the 
progeny of the faint-SGB and suggested that intermediate-mass 
AGB stars might have contributed to the enrichment of the intra 
cluster medium (ICM) before the formation of the second gener- 
ation of stars. Interestingly, both the i-rich and the i-poor groups 
exhibit their own Na-O anticorrelation, which suggests that NGC 
1 8 5 1 has experienced a complex star-formation history. 

lYong & Grundahll (l2008h also suggested the presence of a 
shght metallicity spread (<0.1 dex) among NG C 1851 RGB 
stars. This result has been recently confirmed by ICarretta et al.l 
(120101) on the basis of a larger sample of stars. Fol lowing a clas- 
sificati on scheme based on Fe and Ba abundance. ICarretta et al.l 
(I20T0I) distinguished between a metal-rich, barium-rich (MR) 
and a metal-poor, barium-poor population (MP). They associ- 
ated the MR and the MP components to the bright- and the faint- 
SGB respectively, which is at odds with what was suggested by 
lYong & Grundahil d2008h . 

While the above listed spectroscopic stu dies targeted 
evolved stars in NGC 1851, before this work, onlv lPancino et al.l 
d2010l) analyzed a sample of unevolved stars in this clusters to 
provide index measurements, and we present for the first time an 
abundance analysis of MS, TO and SGB stars. In particular, we 
observed stars located on the faint- and bright-SGB and derived 
for them C and N abundances; aiming for the first time to 
provide insights on the chemical signature diff'erences between 
the two discrete sequences on the SGB. The paper is organized 
as follows: in Sect. |2| we discuss the observations and data 
reductions; in Sect.|3]we define the index passbands and present 
our results from index measurements in Sect. 13.21 Section |4| 
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Fig. 1. Area covered by the catalogs from which the spectro- 
scopic selection was made: the FORS2 in the outer part (shown 
in black) and the ACS/HST field (continuous red line) in the in- 
ner region. The red crosses represent the IMACS+FORS2 sam- 
ple. We also indicate the half tidal radius in this figure (dotted 
circle). 



contains a description of model atmosphere parameters and 
abundance derivations; Sect. I4.3l presents the abundance results; 
Sect. l5.3l focuses on C and N abundances of stars strictly located 
on the faint-SGB or bright-SGB. Finally, we summarize our 
findings and draw the conlusions in Sect. |6] 



2. Observations and data reduction 

2.1. Source catalogs and sample selection 

We selected our targets from lit erature photometry: FORS2 
V and / photometry presented by IZoccali et alj d2009l) . in the 
southwest quadrant of the cluster, as well as F6Q6W and F8 14W 
HST/ACS photometry from the GGC treasury program GO- 
10775 for the inner part of the cluster (M08). The area covered 
by the catalogs and the selected spectroscopic targets is shown in 
Fig.[T] We transformed the coordinates using 2MASS as a refer- 
ence astrometric catalog, so the final catalog is on the same rel- 
ative astrometric system. Spectroscopic targets were selected as 
the most isolated stars located around the turn-off and the SGB, 
reaching the RGB base. The resulti ng photometry was calib rated 
using stars in common with the V, / iBellazzini et al.l d200ll) cata- 
log, covering an 8' x8' field centered on the cluster. We also used 
publ icly available Stromgre n u, v,b,y photometnu of N GC 1851 
from lGrundahl et al.l d 19991) and lCalamida et al.l d2007l) . We refer 
the reader to that paper for details about observations and data 
reduction of Stroemgren photometry. 



http : //www .mporzio . astro . it/spress/stroemgren. php. 
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2.2. Observations and spectroscopic reductions 

We acquired low-resolution (R:^1123 and R:^1246 at 3880 
and 4305 A, respectively) spectra of turnoff and SGB stars in 
the globular cluster NGC 1851 with the IMACS multi-object 
spectrograph at the Magellan 1 (Baade) telescope at the Las 
Campanas Observatory in Chile. The adopted instrumental setup 
with the grating GRAT 600-1 covers the nominal spectral range 
between 3650-6750 Awith a dispersion of 0.38A/pix. This spec- 
tral range includes CN (3880 A) and CH (4305 A) molecular 
bands used to derive nitrogen and carbon abundances. However, 
the actual spectral coverage depends on the location of the slit 
on the mask with respect to the spectral dispersion. Because we 
observed faint stars, the total integration time was long, requir- 
ing ten exposures of 1800 sec each. Therefore all observations 
were made with a single-mask setup with 48 slits. We were able 
to extract spectra for 46 targets from our initial target list. 

To these 46 spectra observe d with Magell a n, we added 47 
other MS and SGB spectra from lPancino et all (1201 Ol) observed 
with the FORS2 multi-object spectrograph on the ESO VLT at 
the Paranal Observatory in Chile. We refer the reader to that 
paper for details of the FORS2 observat ions. We reduced ou r 
data following the procedure described in lPancino et al 
For the data pre-reduction, we used the standard procedure for 
overscan correction and bias-subtraction with the routines avail- 
able in the noao.imred.ccdred package in IRAlQ. Cosmic rays 
were removed with the IRAF Laplacian edge-detection routine 
(Ivan Dokkumll200lh . The frames were then flat-fielded and re- 
duced to one dimension spectra with the task apall. Once we 
obtained ten wavelength-calibrated, one-dimensional spectra for 
each star, we co-added them on a star-by-star basis to reach a 
relevant S/N (typically between 20-30 per pixel at 3880A) even 
in the bluer part of the spectrum 0. As a final step, we examined 
each spectrum and rejected t hose spectra w i th bad quality, fol- 
lowing the recipes outlined in lPancino et al. l (l2OT0h . We defined 
several criteria for this rejection: 

1. S/N ratio < 10 (per pixel) in the CN 3880A band; 

2. clear defects (like spikes or holes) from an individual inspec- 
tion of the spectrum on the band or continuum windows; 

3. discrepant Ca(H-HK) and Hp index measurements. 

Forty-three stars survived our criteria of selection. Moreover, 20 
stars had spectral or continuum passbands falling in the gap be- 
tween the CCDs because of the location of the slit on the mask 
with respect to the dispersion direction. We were not able to 
measure CH and CN band strengths for these stars, but in some 
cases we determined N and C abundances (Sect. |4]i. Therefore 
we decided to retain these spectra and consider them in the sub- 
sequent analysis. 

2.3. Membersiiip and quality control 

NGC 1851 (I,b = 244.51, -35.08;|Hams|[l99l 2010 edition) 
is projected at a far distance from the Galactic plane and conse- 
quently the contamination from field stars is almost negligible. 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

' Before co-adding we checked that the shifts between spectra from 
different exposures were negligible compared to our wavelength cali- 
bration uncertainty (30 km s"') . 
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Fig. 2. Color-magnitude diagrams for NGC 1851. Gray dots 
show V-I photometry for the outer field (left panel) and inner 
field [right panel) from FORS2 and ACS, respe ctively. White 
dots mark spectroscopic targets presented in iPancino et al.l 
while black dots show our newly observed stars. 



In addition, the average radial velocity of cluster stars signifi- 
cantly differs with respect to the field (320.5 km/s; lHarrislll996t 
2010 edition), hence non cluster members can be easily identi- 
fied from their radial velocities. For the stars in the ACS/WFC 
field of view, member stars were furthermore selected on the ba- 
sis of their proper motions (see M08 for details). 

For all the spectra, radial velocities were measured with the 
IRAF task fxcor, which performs the Fourier cross-correlation 
between the object spectrum and a template spectrum (the lat- 
ter with known radial velocity). As a template we chose the 
spectrum of a star with the highest S/N: its radial velocity was 
computed using the laboratory positions of several strong lines 
(e. g. Hp, Ha, Hy and CaHK among others) with the IRAF 
task rvidlines. To derive a robust determination for the radial 
velocities, we performed for a given star the cross-correlation 
against the template in four different spectral regions that span 
the entire spectral coverage, from the bluest part out to the red- 
dest part of the spectrum. Then, the four values were averaged 
together, obtaining typically internal errors of ~ 25-30 kms'K 
We obtained an average Vrad of 317 kms^^ with a dispersion 
of II k ms^^, which fully agrees with the previoii s deter mina- 
tion by Yong & Grundahlj (120081) . IVillanova et al] (|20IO|) . and 
.Carretta et al., (,2O10l) ! 

2.4. Additional literature data 

To increase our observed sample, we added 47 MS and SGB 
stars o bserved with FORS2 that were presented bv lPancino et al.l 
(I2OIOI) . The typical resolution of these spectra is R=/l/5/l^800 
and the S/N ratio (per pixel) in the CN band region is 
S/N(CN)^30. More details concerning th e obser vations and data 
reduction can be found in iPancino et al. and references 

therein. We reduced our data following the same procedure as 
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adopted by these authors. Moreover, their definition of in d ices is 
exactly the same as ours (Sect. O. While IPancino et aH (1201 Ol) 
did not need to normalize their spectra, this was necessary in 
our case, so we fitted a polynomial to the pseudo continuum 
to normalize both IMACS and FORS2 spectra as described in 
Sect. [3] We adopted the same procedure in determining the C 
and N abundance for both the newly obtained and previously 
studied spectra. 



3. Index definition and measurements 

Spectral indices are defined as a window centered on the molec- 
ular band of interest and one or two windows around it to define 
the continuum level. For each spectrum, S3839 and CH4300 in- 
dices sensitive to the absorption by the 3883A CN band and the 
4300A CH band were measured. Several spectral index defini- 
tions exist in literature, generally optimized to quantify the CN 
content of the atmospheres of red gia nt stars. In our case, to be 
consistent with the previous work of iPancino et al. (2 010|). w e 
decided to adopt the indices as defined in lHarbeck et al. 



53839(CA^) = -2.5 log 



-^3861- 



3884 



54142(CA^) = -2.5 log 



Ci/4300 = -2.5 log 



■^3894-3910 
■^4120-4216 



0.5/^4055-4080 + 0.5^4240-4280 

■^4285-43 15 



0.5/^4240-4280 + O.5F4390-4 



In part i cular, the S3839(CN) index we use d differs from that of 
iNorrislJTgsTh oriN orris &Freemanl (119791) . and it accounts for 
stronger hydrogen lines in the region of CN feature for stars 
cooler than red giants. As an additional check, we defined and 
measured two different CN band indices in the wavelength re- 
gion covered by our spectra: the S3839 for the CN band around 
3880 A and the S4142 for the one around 4200 AQ. We ob- 
tained in dex measurement uncer t ainties with the expression de- 
rived by IVollmann & Eversberd (l2006l) . assuming pure photon 
noise statistics in the flux measurements. In addition, we mea- 
sured the two indices centered aro und the calcium H an d K lines 
and the Hp line as defined again in IPancino et al ] (l2010h to reject 
outliers from our sample. 

The final reduced spectra generally show a strong decline 
of the signal toward bluer wavelengths. This is largely expected 
and may be due to the different instrumental efficiency, higher 
absorption of the Earth's atm osphere in the blue and stellar flux 
wavelength dependency (see ICohen et al.ll2002l for a complete 
discussion). The CH band at 4300 A is not affected by the 
change in spectral slope from atmosphere or instrumental effects 
thanks to two continuum bandpasses. On the other hand, we had 
to rely only on a single continuum bandpass in the red part of the 
spectral feature for the 3883A CN band. Following lCohen et al.l 
(2002, 2005), we decided to normalize the stellar continuum in 
the spectrum of each star, then found the absorption within the 
CN bandpass. Moreover, by fitting the continuum, we were able 
to directly compare the indices measured in this section and the 
abundances derived from spectral synthesis in Sect|4] We fitted 



* As discussed by different previous authors, the S3839 index is 
found to be by far much more sensitive to CN variations with respect 
to S3839. For the S4142 index, the spread is generally of the size of (or 
slightly wider) than the median error bar on the index measurements. 
Therefore, we decided to rely only on the S3839 index measurements 
throughout. 




4200 

wavelength (A) 




4300 4320 4340 

wavelength [k] 



3850 3900 3950 

wavelength [A] 



Fig. 3. Panel a: example of the windows from which we mea- 
sured the CN and CH indices (dark gray hatched regions) to- 
gether with their respective continuum windows (light gray 
hatched regions). Panels b and c show the windows adopted 
for the H and K calcium index and the Hp index. The non- 
normalized superimposed spectrum (star 41213, S/N~35 in the 
S3839 region) was smoothed for clarity. Panel d: Observed 
(small black dots) and synthetic (line) spectra around CH band 
for the star 15490. The black best fits, while the red and green 
lines are the syntheses computed with C abundance altered by + 
0.10 dex from the best value. Panel e: The same as in panel d 
but for the CN feature for the star 40062. The synthetic spectra 
show the best fit (thick black line) and the syntheses computed 
with N abundance altered by + 0.20 (thick red and green lines). 



a thi rd-order polynomi al masking out the region of the CN band 
(see ICohen et al.ll2005h . The polynomial fitting used a 60" high 
and 3(T low clipping, running over a five pixel average. Then 
we computed S3839 indices from these continuum-normalized 
spectra and used the average (0.126+0.04 and 0.05+0.01 for 
IMACS and FORS2 spectra respectively) to set a zero point off- 
set and thus delete the instrumental signature present in the raw 
S3839 indices. 

The measured indices are listed in Tab. [1] and plotted in 

Fig. a 



3.1. Dependency on temperature and gravity 

CN and CH bands are stronger at a fixed overall abundance 
for stars with lower temperature and gravity. In particular, the 
formation efficiency of the CN molecule strongly depends on 
the temperature, therefore we expect the indices depend on the 
color of t he single stars. These dependencies are usually elimi- 
nated (see lHarbeck et al.ll2003HKavser et al.ll2008L for example) 
by fitting the lower envelope of the distribution in the index- 
magnitude plane (or index-color plane). For our sample, we used 
the median ridge line, shown as d ashed red lines in F ig. H] to 
eliminate these dependencies (see IPancino et al. l l2oToh . These 
rectified CN and CH indices are in the following indicated as 
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Fig. 4. Measured S3839 and CH4300 indices fo r the program 
stars a s a function of V. White dots are stars from lPancino et all 
(120 lOh . while black dots are from the IMACS sample. Dashed 
red lines in both panels are the median ridge lines used to ehmi- 
nate temperature and gravity dependencies. 



Fig. 5. Plot for the distribution of CN and CH band strengths. 
Gray dots show measurements for stars. CH-weak and CH- 
strong stars are separated by the horizontal dashed line and their 
centroids with error bars (drawn at lo") are marked as black dots 
and large empty dots, respectively. 



5S3839 and 5CH4300, respectively, and we refer to these new 
indices throughouj^. For the stars in common between this work 
and Pancino et al. (2010), the mean difference in the 5S3839 and 
(JCH4300 indices, derived by subtracting Pancino et al.'s values 
from ours, are 0.00 + 0.02 mag and -0.01+0.01 mag, respec- 
tively. Because we a dopted the same reduction procedures as in 
iPancino et al.l (l2010l) . we can only ascribe this small difference 
to the continuum normalization we performed (see Sect.[3]l. 



3.2. CN and CH distribution 

Variations of several light elements and anticorrelations be- 
tween strengths of the CN and CH bands were detected for very 
many clusters. Because molecular abundance is controlled by 
the abundance of the minority species, the corrected CH index 
is a proxy for carbon abundance, while 5S3839 traces the nitro- 
gen abundance. The visual inspection of the top panel of Fig.|4] 
reveals significant scatter in the CN index over the magnitude 
range with V > 19.5 with some hints of bimodality toward the 
brightest tail of the distribution. The range of CN becomes less 
evident at fainter luminosities as a consequence of the increas- 
ing temperature. In the bottom panel of the same figure, we plot 
the CH4300 and S3839 versus the stellar V-band magnitudes. 
Here the variations among the measured index are very small 
and within the uncertainties throughout. 

Figure |5] shows the rectified index 5S3839 as a function 
of 5CH4300 for all stars. We found no evidence for a signif- 



^ We obtained a rough estimate of the uncertainty in the placement of 
these median ridge lines by using the first interquartile of the rectified 
indices divided by the square root of the total points. The uncertainties 
(typically ~ 0.01 for the CN index and ~ 0.005 for the CH index) are 
largely negligible for the applications of this work. 



ica nt CH-CN antico rrelation, similarly to what was found by 
Pancino et all (1201 Ol) . 



4. Spectral synthesis and abundance derivations 

Indices are a fast tool to characterize chemical anomalies, but 
we can also rely on spectral synthesis to fully characterize our 
target stars. This becomes necessary when indices do not offer 
conclusive answers, as we saw in the previous sections. 



4.1. Atmosplieric parameters 

We derived estimates of the atmospheric parameters from the 
calibrated ACS and FORS2 photometry presented in Sect. |2] 
Dere ddened ( V - /)o colors were obtained adopting E{B - V) = 
0.02 (lHarriJl996l 2010 edition). We obtained effective temper- 
atures and bolometric corrections (hereafter Tgff and BCy) with 
the Alonso et al. ( 1996, 1999, 2001 ) colo r -temp erature relations, 
adopting [Fe/H]= -1.22 from lYong et^ (l2009l) . and taking into 
accou nt the uncertainties in the m agnitudes and reddening esti- 
mates. lAlonso et al.l (1 19961 Il999h adopted Johnson's system as 
a reference, therefore we converted {V - Ic) into (V - Ij) after 
der eddening using the prescriptions by |Bessell (197^ to feed 
the lAlonso et al.1 (Il996[ll999l 1200 Ih calibration. Gravities were 
then obtained by means of the fundamental relations 

log^=log— +21og^, 
0.4 {Mhoi - MhoiQ) = -4 log + 2 log ^® 



R 



wh ere we assumed the solar values reported in lAnderseiil 
(119991) : log^o = 4.437, T.ff^ = 5770K and Mtoi,Q = 4.75. 
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For all our stars, we assumed a typical mass of 0.8 Mq 
(iBergbusch & Van denBerg' '2001*) and a distance modulus of 
(m - M)v = 15.47 (Harris 1996, 2010 edition). 

Finally, we obtained the microturbulent velocities (v,) from 
the relation log^ a nd v,, i.e., v, - 1.5 - 0.03 logg as in 
ICarretta et al.l (l2004l) . This method leads to an average micro- 
turbulent velocity estimate of v, = 1.0 + 0.1 km.?"', therefore 
we assumed v, = 1.0 km i"' for the entire sample. An addi- 
tional check to test the reliability of our atmospheric param- 
eter determination was performed usin g theoretical isochrone s 
downloaded from the BaSXfl database (IPietrinferni et al ] l2006h . 
We chose an isochrone of 11 Gyr (12 Gyr for the faint-SGB) 
with standard a-enhanced composition, and metallicity Z = 
0.002 and we projected our targets on the isochrone to obtain 
their parameters. We present the average difference between 
the two methods in Fig. |6] In the top panel we plo t the d iffer- 
ence in temperature obtained using the I Alonso et al ] (119991) cal- 
ibration and the temperature obtained by projection of the star 
on the BaSTI isochro nes as a funct i on of the temperature de- 
rived by means of the lAlonso et al.l (Il999h empirical relations. 
The scatter for high temperatures is significant, as largely ex- 
pected, because these empirical calibrations were obtained for 
giant stars and therefore are valid in a precise range of color. On 
the other hand, the scatter is modest (and in many cases within 
uncertainties) when comparing temperatures obtained by using 
th^lonso et al. (1999) and Alonso et al. (1996) relations (see 
bottom panel of Fig. |6]l, where the latter was derived for low 
main sequence stars. However, we preferred to avoid using tem- 
peratures derived by isochrone fitting mainly for these reasons: 
(a) we cannot assume a priori that the cluster is a single popu- 
lation (with the same [Fe/H] and CNO content, among others) 
and (b) the projection on the {V, V-I) plane is always uncertain, 
and a rigorous treatment should include (asymmetrical) errors 
on the V magnitude and V-I color, and finally, (c) different 
sets of isochrones (Padova, BaSTI, and DSEP for example) give 
different results. 

As discussed in the following sections, even if the differences 
between the two temperatures scales appear non-negligible, the 
main results of the paper appear totally unchanged if we adopt 
one or the other temperature scale. This is mainly because the 
abundances ranking among targ et stars is left unch anged. We 
therefore preferred to rely on the lAlonso et al. parameter 
estimates and discuss the effect of the chosen temperature scale 
below. 



4.2. Abundances derivation 

We use d the local the rmodynamic equilibrium (LTE) program 
MOOG (ISnedenI 19731) combi ned with the ATLAS9 model atmo- 
spheres (iKurucdl 19931 l2005h to determine carbon and nitrogen 
abundances. The atomic and molecular line lists were taken from 
the latest Kurucz compilation and downloaded from F. Castelli's 
websiteQ. 

Model atmospheres were calculated with the ATLAS9 code 
starting from the grid of models available in F. Castelli's web- 
site, using the values of Tg/f, \ogg, and v, determined as ex- 
plained in the previous section. The ATLAS9 models employed 
were computed with the new set of opacity distribution func- 
tions (ICastelli & Kuruczl2003l) and excluding approximate over- 
shooting in calculating the convective flux. For the CH transi- 
tions, the log g obtained from the Kurucz database were revised 
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Fig. 6. Top panel: Differences in the temperature estimates 
by means of isoch rone projection procedure and using the 
Alonso et al. (1999) relation (A(r,,o - TAg9)) as a function of the 
Alonso et al. ( 1999J temperature (7^ 99) for all our target stars. 
Bottom panel: The same, but for the lAlonso et alj d 19961) caU- 
bration (r^96). 



do wnward by 0.3 dex to b etter reproduce the solar- flux spectrum 
bvlNeckel & Labsl (11984 with the C abundance bvlCaffau et al.l 
( I2OII ). as extensively discussed in lMucciarelli et aiT ( 201 lb . 

C and N abundances were estimated by spectral synthesis 
of the 2S-2n band of CH (the G band) at -431 OA and the UV 
CN band at 3883A (including a number of CN features in the 
wavelength range of 3876-3890A), respectively. Lower panels 
of Fig. [3] illustrate the fit of synthetic spectra to the observed 
ones in CH and CN spectral regions. Abundances for C and N 
were determined together in an iteractive way, because for the 
temperature of our stars, carbon and nitrogen form molecules 
and as a consequence their abundances are related to each other. 
The input model atmosphere was used within the MOOG run- 
ning synth driver that computes a set of trial synthetic spectra 
at higher resolution (0.3A intervals) in the spectral region be- 
tween 41 50-4450A, varying the carbon abundance in steps of 0. 1 
dex typically in the range of -0.2 to -1.2 dex to fit a full band 
profile. After the synthesis computations, the generated spectra 
were convolved with Gaussians of appropriate FWHM to match 
the resolution of the observed spectra. In this way, the carbon 
abundances were derived by minimizing the observed-computed 
spectrum difference and were used to determine A(C). The car- 
bon abundance was then used as input in the synthesis of the 
UV CN feature to derive nitrogen abundances. The procedure 
was repeated until we obtained convergence within a tolerance 
of 0.1 dex in the C and N abundances. 

For the results presented here, a fit was determined by 
minimizing the observed-computed spectrum difference in a 
60A window centered on 4300A for the CH G-band and 
40A window for the UV CN feature at 3883A. Running synth 
on quite a broad spectral range (200 and 300A for the G-band 



http : //wwwuser . oat . ts . astro . it/castelli/linelists . html^ d the CN feature, respectively) to produce synthetic spectra 
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Fig. 7. [C/H] and [N/H] abundances for the NGC1851 SGB, TO 
and MS stars in Table |2] are plotted. A C versus N anticorrela- 
tion is evident. Stars that have already experienced some mixing 
episode are shown as large white dots, while three notable out- 
liers are represented as large triangles (see text for comments). 
Histograms of [N/H] and [C/H] with typical median error bars 
are also plotted in this figure. 

allowed us to set a reasonable continuum level also by visual 
inspection and thus compute robust abundances. 

We adopted a constant oxygen abundance ([O/Fe]=0.4 dex) 
throughout all computations. The derived C abundance is de- 
pendent on the O abundance and therefore so is the N abun- 
dance, and in molecular equilibrium an over-estimate in oxy- 
gen produces an over-estimate of carbon (and vice versa), and 
an over-estimate of carbon from CN features is refleced in an 
under-estimate of nitrogen. We expect that the exact O values 
will affect the derived C abundances only negligibly, since the 
CO coupling is marginal for stars warmer than ~ 4500 K. To 
test the sensitivity of the C abundance to the adopted O abun- 
dance we varied the oxygen abundances and repeated the spec- 
trum synthesis to determine the exact dependence for a few rep- 
resentative stars in a wide range of T^fj (from 5200 to 5900 K). 
In these computations, we adopted [0/Fe]= -0.5 dex, [0/Fe] = 
0.0 dex, [0/Fe]= +0.5 dex. We found that strong variations in 
the oxygen abundance slightly affect {5A{C) 1 5{0 1 F e\ ^ 0.15 
dex) the derived C abundance in colder stars iTeff< 5400 K), 
while they are completely negligible (on the order of 0.05 dex or 
less) for warmer stars. This is within the uncertainty assigned to 
our measurement. 

The total error in the A(C) and A(N) abundance was com- 
puted by taking into account the two main sources of uncer- 
tainty: (i) the error in the adopted T^ff, typically dA{C)/dTeff - 
0.08-0.10 dex and 6A{N)/6Teff - 0.11-0.13 per 100 K for the 
warmest stars in our sampl^l; (ii) the error in the fitting proce- 
dure and errors in the abundances that are likely caused by noise 

Cooler stars are slightly less sensitive to T^ff variations, typically 
at a level of6A{C)/6T,ff ^ 0.07-0.08 dex and 6AiC)/6T,ff ^ 0.09-0.1 1 
dex per 100 K. 



in the spectr^. The errors due to uncertainties on gravity and mi- 
croturbulent velocity are negligible (on the order of 0.02 dex or 
less). The sensitivity of the derived abundances to the adopted at- 
mospheric parameters was obtained by repeating our abundance 
analysis and changing only one parameter at each iteration for 
several stars that are representative of the temperature and grav- 
ity range explored. Thus, we assigned the internal error to each 
star depending on its T^ff and \ogg. The errors derived from the 
fitting procedure were then added in quadrature to the errors in- 
troduced by atmospheric parameters, resulting in an overall error 
of ~ +0. 14 dex for the C abundances and ~ +0.28 dex for the N 

values. 

Very recentlv ICarretta et alJ (1201 Ol) found in NGC 1851 a 
small spread in metallicity for a large number of giants, which 
is compatible with the presence of two different groups of stars 
whose metallici ty differs by 0.06-0.0 8 dex (but this result was 
not confirmed in lVillanova et al.l2010l) . This finding could affect 
our analysis in principle, because we adopted the same metal- 
licity for all our stars in the synthesis ([Fe/H]=-1.22). To test 
this effect we repeated the synthesis by altering the metallic- 
ity of stars belonging to the faintest SGB by 0. 10 dex (well 
above the spread claimed bv ICarretta et aljr2010l) . The result- 
ing abundance variations are within the uncertainty assigned 
to our measurement (typically 6A{C)/d[Fe/H] ^ 0.07 dex and 
6A{N) 1 6{F e I H] ^ 0.04 dex) for our low-resolution spectra. 
Therefore this potential small [Fe/H] variation among our spec- 
troscopic targets has no influence on our analysis and conclu- 
sions. We present the abundances derived as described above 
and the relative uncertainties in the abundance determination in 
Tab. |2] Additionally, this table lists the derived atmospheric pa- 
rameters of all our targets. 

4.3. C and N abundance results 

Fig. Qplots [C/H] versus [N/H] measured in this paper. We ob- 
serve strong star-to-star variations in both elements, as already 
observed in all GCs studied to date. An anticorrelation, with con- 
siderable scatter, is apparent from Fig. |2l The scatter is consis- 
tent with the observational errors, but there are a few outliers. In 
a sample of 64 objects with Gaussian errors, two outliers at the 3 
cr level are not expected. The deviation of stars 41350 (V = 19.3) 
and 40022 (V=19.9) with extremely depleted C, from the mean 
relation shown by the NGC 1851 sample in FigQis of higher 
statistical significance. We cannot provide a reliable exp lanation 
for this. Both stars are from the iPancino et al.l (1201 Ol) sample 
and, judging from their radial velocities, are cluster members. 
Moreover, their V,I magnitudes do not have large errors. As a 
tentative hypothesis we suggest that these two stars could be- 
long to the extreme population, using the scheme suggested by 
ICarretta et all ( l2009h . 

All our stars are C-depleted, with moderately weak varia- 
tions in carbon abundances (from [C/H]~-2.7 to [C/H] — 1.5) 
antic orrelated with strong variations in N. The nitrogen abun- 
dance spans almost 2 dex, from [N/H] — 1.9 up to [N/H]~0.0 
de?jPi To check the dependence of the carbon and nitrogen abun- 
dances on the adopted temperature scale, we re-ran the synthesis 
using the atmospheric parameters derived by isochrone fitting 



' Additionally, in the treatment of internal error for nitrogen we var- 
ied the carbon abundance by ±0.10 dex (that is the typical error asso- 
ciated to A(C)). We added these errors in quadrature with those intro- 
duced by the model atmosphere to estimate the internal uncertainty of 
the A(N) values. 

if we exclude the outliers discussed above. 
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Fig. 8. Comparison of carbon (top) and nitrogen (bottom panel) 
abundances derived by adopting differ ent temperature scales 
("A99" refers to the lAlonso et"al1l 19991 calibration and "iso" to 
the isochrones fitting procedures) with their relative uncertain- 
ties. 



(see Sect. 14.11 ). The result of this exercise is shown in Fig. |8] 
As can be seen from this figure, the carbon abundances would 
be higher considering these higher temperatures (ranging from 
~ 0.2 dex for giants up to 0.4-0.5 dex for Ms stars). This reflects 
on the nitrogen abundances, as demonstrated in the bottom panel 
of the same figure. This is only to show that the abundances rank- 
ing among target stars is left unchanged: while the zero point of 
our derived abundances would shift, the amplitude of the star- 
to-star variations for C and N would remain similar regardless 
of the adopted stellar parameter Therefore, our conclusions do 
not depend on the adopted stellar parameters. 

Evaluating the accuracy of our absolute abundance scale is 
very difficult because we found no literature data to compare 
with. Fig. |9] compares the C and N abundances of this paper to 
the abundances derived for M 5 bv lCohen et all (l2002h. a clus- 



ter with a metallicity comparable ( — 1.29: lHarrislll996l 2010 
edition) to that of NGC 1851. In panel (a) we present a com- 
parison between C and N abundances derived by assuming the 
lAlonso et al.l (1 19991) te mperature scale and abundances derived 
bv lCohen et al.l (l2002h for stars at the base of the RGB in M 5, 
while panel (b) refers to isochrones-fitting temperatures. In both 
panel a clear C-N anticorrelation is apparent. According to theo- 
retical computations and earlier investigations, the carbon abun- 
dance declines fro m MS to RGB . In pan el (a) there is a mild dis- 
agreement with the ICohen et al.l (l2002l) data; that is completely 
recon ciled in pane l fZ?). T his effect can be entirely explained be- 
cause jCohenetal] (|2002) used atmosphere parameters obtained 
from the isochrone. Still, although there is an offset, the two 
anticorrelations seem to follow a similar pattern. We conclude 
again from Fig. |9|that the anticorrelation we observe is totally 
untouched by the choice of the temperature scale, and shifts in 
the absolute abundance scale cannot account for the wide range 
in N abundances apparent in FigjT] 



Fig. 9. Panel (a): [C/Fe] and [N/Fe] abundances for NGC 1851 
stars. Stellar atmospheres and spec tral syntheses were derived by 
assuming the lAlonso et al.l (Il999l) temperature scale. The sym- 
bols are the same as in Fig. U] Abundances and relative uncer- 
tainties for stars in M 5 from lCohen et a are also shown 
as crosses for comparison. Panel (b): the same as in the left 
panel, but assuming temperatures and gravities obtained from 
the isochrones. 



We therefore conclude that the C versus N anticorrelation 
among unevolved NGC 1851 stars in FigjT] is indeed real and 
from here on we wi ll therefore only present results based on the 
lAlonso et al.l d 19991) temperature scale. 

We also plotted the derived abundances as a function of the 
V magnitude and V - I color in Fig.[TO|to evaluate possible sys- 
tematic effects with luminosity and temperature. 

While none of these effects are apparent, we can tentatively 
identify the occurrence of a mixing episode for NGC 1851 stars 
from this plot. The top panel of Fig. [TOl shows a notable decline 
in the carbon abundances for stars with V <18.9 and (V-I) >0.8 
(stars marked as white dots in the same figure), which is ex- 
pected for stars in the course of normal stellar evolution. This 
behavior of the C abundance allows us to identify stars that ex- 
perienced a major mixing episode, which may alter the primor- 
dial abundances. Curiously enough, these stars, plotted again as 
large white dots, seem to define a pretty clear and narrow anticor- 
relation in Fig|7](Spearman's rank correlation coefficient -0.92). 
The shape of this anticorrelation agrees with what we expect af- 
ter the occurrence of a mixing episode: the high N enhancement 
found in unevolved or less-evolved stars is strongly softened by 
evolutionary effects and a large part of dwarfs and early sub- 
giants have N abundances as high as those observed in slightly 
evolved RGB stars. We identify, again from Fig.[TO| three out- 
liers, coded as empty triangles as in Fig. [Tj Two of these stars 
were found to deviate significantly from the main C-N relation 
if Fig. I2I We call these anomalous only by virtue of their posi- 
tions in the upper panel of Fig. [TOjand decided to not consider 
them further 
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Fig. 10. Derived C and N abundances are plotted against the pho- 
tometry for NGC 1851 stars. No systematic trends with either lu- 
minosity or temperature are apparent in the abundances. We call 
stars marked as black empty triangles anomalous, while stars 
that supposedly underwent some mixing episodes are plotted as 
empty dots (see text). 



At this point we note that we cannot arbitrarily distin- 
guish between two groups of stars with different [N/H] or 
[C/H] because we are unable to detect any clear bimodality. 
To be more quantitative, w e ran the dip test on unimodality 
dHartigan & Hartiganlll985b . We performed this simple statisti- 
cal test only on stars with a magnitude V <18.9 and can confirm 
that there is no bimodality in either the [C/H] or [N/H] . 



5. The chemical composition of the double RGB 
and SGB 

As already discussed in Sect. [T] the discovery of multiple se- 
quences in the CMD of NGC 1851 provided unambiguous proof 
of the presence of multiple populations and brought new interest 
and excitement about this GC. While it is now widely accepted 
that NGC 1851 hosts two or more stellar populations, the con- 
nection among its multiple SGBs, RGBs, and HBs is still con- 
troversial and the chemical composition of the two SGBs is also 
debated. 

Several authors suggested that the groups of 5-rich and 
s-poor stars detected from RGB-star spectroscopy are the 
progeny of the faint- an d bright-SGB, respectively (e. g. M08, 
lYong & Grundahll l2008h . in close analogy to what was ob- 
served in M22 and co Centau r i (e. g. iMarino et al.ll2009l 120111: 
Johnson & Pilachowskil l20l"ot iPancino et al.l 201 1 ). In contrast. 
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Carretta et al. (2011) claimed that the faint-SGB consists of 
barium-poor metal-poor stars while bright-SGB stars have an 
enhanced barium and iron abundance. 



1 2 3 4 5 6 

(>j -y) + (v-h) 

Fig. 11. CMD for NGC 1851 from Stromgren photometry. The 
inset shows a zoom of the SGB region. Only stars with high- 
quality photometry were plotted (see Sect. l5.1l for details). Note 
the discrete double RGBs connected with the bimodal SGBs. 



5.1. Photometric connection between SGB and RGB 

To investigate this question in more detail, we started an- 
alyzing literature photometry. We used the WFC/ACS HST 
CMD in F606W and F814W bands pres ented in M08 (see 
ISaraiedini et al.ll2007l: lAnderson eTai1l2008i for deta i ls) and the 
Stromgren u, b, v, y ph otometry from lGrundahl et al.l ( 1 19991) and 
ICalamida et alJ (l2007h . Here we are interested in high-quaUty 
photometry and included in the analysis only relatively isolated, 
unsaturated stars with good values of the PSF-quality fits and 
small rms errors in astrometry and photometry. A detailed de- 
scripti on of the selection procedures is given in iMilone et al.l 
(I2OO9I Sect. 2.1). We corrected our photometry for remaining 
spatially depen dent errors, c aused by small inaccuracies of the 
PSF model (see lAnderson et al. 2008). To account for the color 
differences of these variations we followed the recipes from 
iMilone et al.l (1201 ll Sect. 3). Briefly, we defined a fiducial line 
for the MS by computing a spline through the median colors 
found in successive short intervals of magnitude, and we iter- 
ated this step with a sigma clipping; then we examined the color 
residuals relative to the fiducial and estimate for each star, how 
the observed stars in its vicinity may systematically lie to the 
red or the blue of the fiducial sequence. Finally we corrected the 
star's color by the difference between its color residuals. 

To study multiple populations from the CMD analysis, we 
started searching for the combination of magnitude and colors 
that provides the best separation of the two RGBs and SGBs in 
NGC 1851. Results are illustrated in Fig. [TTI where we plot u as 
a function of (m -y)+{v - b). A visual inspection of this diagram 
leaves no doubts on the presence of a bimodal RGB and SGB and 
shows that the faint-SGB and the bright-SGB are clearly con- 
nected with the red- and the blue-RGB, respectively. A similar 
connection between the two SGBs and RGBs has already been 
observed for NGC 1851 bv iHan et all (l2009l) i n the U versus 
(U — I) CMD and was studied more recently bv lSbordone et al.l 
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Fig. 12. Main panel: Reproduction of the CMD of Fig. [TT] 
We used red color codes for red-RGB and faint-SGB stars, 
while blue-RGB and bright-SGB stars are represented in blue. 
Symbols refer to stars from spectroscopic studies as indicated in 
the figure. The rectified CMD and the histogram color distribu- 
tion of RGB stars between the two dashed lines is shown in the 
inset (see text for details). 

(1201 Ih . These authors showed that while the double SGB is con- 
sistent with two groups of stars with either an age difference of 
about one Gyr or with different C+N+0 overall abundance, the 
double RGB seems to rule out the possibility of a large age dif- 
ference. 

Fig. [TT] revealed that the bimodality found in the SGB can 
also be seen in the RGB. To further confirm our association of 
the faint-SGB (bright-SGB) component with the red-RGB (blue- 
RGB), we focused on the relative number of stars of all evolu- 
tionary stages. To estimate the fraction of stars in the two RGBs 
we used only RGB stars between the two dashed lines in the 
magnitude interval where the split is more evident (main panel 
of Fig. [T2I 1. The procedure is illustrated in the inset of Fig. [121 
To obtain the straightened RGB of the right-hand panel, we sub- 
tracted from the color of each star the color of the fiducial se- 
quence at the M magnitude of the star The color distribution of 
the points in the middle panel were analyzed in two magnitude 
bins. The distributions have two clear peaks, which we fitted 
with two Gaussians (red for the red-RGB and blue for the blue- 
RGB). From the areas below the Gaussians, 70+3% of stars turn 
out to belong to the blue-RGB, and 30+3% to the red one. With 
the statistical uncertainties these fractions are the same in both 
magnitude intervals and roughly match the relative frequency on 
the fainter/brighter SGBs (35% versus 65%, M08) and HB stars 
on the blue/red side of the instability strip (35% versus 65%). 

5.2. Chemical composition of NGC 1851 subpopulations 

Because the chemical abundance determinations presented so far 
do not define any clear bimodality, the clear separation of the se- 
quences of Fig.[TT]provides a unique opportunity to obtain infor- 



Fig. 13. Upper panel: Na-O anticorrelation among NGC 1851 
RGB stars from HR spectroscopy studies. Red color refers to 
stars photometrically selected to belong to the red-RGB in the 
Stromgren m, (m -y) + (v - b) diagram, while stars located on the 
blue-RGB are shown in blue. Symbols and color code are con- 
sistent with those of Fig. [12] Bottom panel: The run of [La/Fe] 
versus [Fe/H] and the normalized number distribution for red 
and blue stars in this plane. The color code is consistent with the 
upper panel. 

mation on the chemical differences between the two RGBs and 
SGBs in NGC 1 85 1 . To do this, we used a m, (m - y) + ( v - fe) di- 
agram to isolate the samples of blue-RGB and bright-SGB stars, 
and red-RGB and faint-SGB stars. Then we plotted with red and 
blue symbols the red-RGB and blue-RGB stars for which abun- 
dance measurements are available from high-resolution spec- 
troscopy. 

Our analysis of the chemical abundance patterns of the two 
RGBs is summarized in Fig. [13] Lower panels show [Fe/H] ver- 
sus the abundances of the ^-process ele ments barium and lan- 
thanum mea sured by lYon g & Grundah| ( |2008|) . IVillanova et al.l 
(1201 Oh . and [Carretta etalj' (l201ll) from GIRAFFE and UVES 
data. The histogram of the s-element distribution is illustrated in 
the middle panel, while upper panels plot [Na/Fe] versus [O/Fe]. 
The average iron, barium, lanthanum, sodium, and oxygen abun- 
dances are listed in Table|3]for the two groups of stars. 

In the light of our analysis of literature photometric and spec- 
troscopic data we are now able to characterize the two RGBs and 
SGBs of NGC 1851 as follows: 

- Faint-SGB and red-RGB stars are photometrically con- 
nected, therefore they represent the same subpopulation of 
NGC 1851; the same can be said abo ut bright-SGB and blue- 
RGB (see also iMarino etani2012bl for the case of M 22). 
This connection is supported by the relative (percentage) 
numbers of the sequ ences; therefore the d ata do not support 
the interpretation bv lCarretta et all (1201 Ol) that the red-RGB 
is associated to the bright-SGB. 

- Literature data suggest that the red-RGB stars tend to be en- 
riched on average in Na and .^-process elements, and poor in 



10 



Lardo et al.: C and N abundances of MS and SGB stars in NGC 1851. 
Table 3. Mean abundances for NGC 1851 stars from high-resolution studies. 



Element 


Abundance (blue-RGB) 




Abundance (red-RGB) 




References 


[La/Fe] 


0.27+0.02 


5 


0.61+0.05 


3 


1 


[Na/Fe] 


-0.05±0.11 


5 


0.57+0.15 


3 


1 


[0/Fe] 


0.50±0.04 


5 


0.17+0.14 


3 


1 


[Fe/H] 


-1.29+0.04 


5 


-1.23+0.07 


3 


1 


[Ba/Fe] 


0.09+0.02 


8 


0.52+0.03 


7 


2 


[Na/Fe] 


0.04+0.11 


8 


0.47+0.07 


7 


2 


[0/Fe] 


0.09+0.07 


8 


-0.19+0.08 


7 


2 


[Fe/H] 


-1.23+0.01 


8 


-1.22+0.01 


7 


2 


[Ba/Fe] 


0.43+0.02 
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References. (l) lYong & Gnindahll l [200g) : (2) IVillanova etal] l [20Tol) : (3) ICarretta et"ai] Jlonl) . 

Notes. Red- and blue-RGB stars are defined according to their location with respect to the ridge line used to define color residuals in Fig. [12] 
GIRAFFE data <*' UVES data 




V-I (u-y) + (v-b) 



Fig. 14. Left panel: selected bright-SGB (marked as blue dots) 
and faint-SGB (shown as red squares) stars are overplotted on 
the y, V-I CMD presented by M08. Right panel: selected bright- 
and faint-SGB stars are shown in the Stromgren m, {u-y) + {v-b) 
diagram. The color code is consistent with the left panel. 
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Fig. 15. Photometrically selected bright-SGB (marked as blue 
dots) and faint-SGB (shown as red squares) stars are represented 
in the A(C) versus A(N) plane. Stars that we were unable to 
unambiguously associate with either of the two populations are 
shown with open circles. 



oxygen, while blue-RGB stars appear to have their own, ex- 
tended anticorrelation and to be solar in Ba and .s--process 
eleme nts. This is particularly evident in the Carretta et al] 
(1201 Ih dataset, which also has the highest statistical value. 

- Red-RGB (and thus faint-SGB, according to our interpre- 
tation above) stars are enhanced in barium and lanthanum 
by ~0.3-0.4 dex with respect to the blue-RGB (and conse- 
quently the bright-SGB). 

- The literature data suggest that there is no significant iron 
difference between the two groups of stars. In this context 



we rec all that lYong & Gnindahll (l2008h and ICarretta et aH 
(1201 Ih detected a significant [Fe/H] variation among both 
i-rich and i-poor stars but these results strictly disagree with 
the narrow iron distribution observed by Villanova and col- 
laborators. The presence of an intrinsic iron spread among 
NGC 1851 stars is still controversial. 
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A(C + N) A(C + N) 



Fig. 16. Left panel: Histogram of the C+N distribution for 
SGB stars selected as described in Sect. 15.31 The median error 
bar is plotted below the histogram for reference. Right panel: 
Histograms of the C+N sum for bright-SGB (dashed blue) stars 
and faint-SGB (hatched red) stars are shown. 

5.3. C and N abundances along the double SGB 

In this section we present the chemical composition of stars on 
the two SGBs of NGC 1851. A bona fide sample of stai's that 
belong unambiguously either to the faint-SGB or to the bright- 
SGB were selected using both the V,V -I and u, (u-y) + (v - b) 
diagrams (see Fig.[T4b. 

The C and N abundances of these bona fide stars are plot- 
ted in Fig. [15] where it is immediately clear that stars belonging 
to the bright-SGB show a fully developed anticorrelation, while 
stars belonging to the faint-SGB appear to have a smaller scatter, 
and to have on average an excess of N (this remains still valid 
when considering temperatures derived by isochrone fitting as 
input of the synthesis, as anticipated in Sect. 14.3b . This new re- 
sult supports our previous identification of the faint-SGB as the 
parent population of the red-RGB, and of the bright-SGB as the 
parent of the blue-RGB, not only on the basis of photometry and 
po pulation ratios , but a lso o n the basis of c hemical composition. 

ICassisi et aP (l2008h and lVentura et al.l Q009 ) argued that the 
overa ll CNO abundanc e difference can account for the SGB 
split. lYong et al] ( l2009l) found evidence for str ong CNO vari- 
ations in contradiction with the the results of IVillanova et alJ 
(l2010l) . To further investigate this hypothesis, we computed the 
C-l-N sum for our bona fide faint-SGB and bright-SGB stars. 
We then derived the histograms of the distribution of the C+N 
sum. In the left panel of Fig. [16] we plot the histograms for the 
entire dataset (faint-SGB + bright-SGB bona fide stars) with 
the typical (median) error bar indicated. The histogram shows 
a high dispersion with a hint of bimodality (with two clumps 
separated at A( C+N)^7.35). For thi s larger dataset, according 
to a KMM test d Ashman et al.l 1 1994) . a bimodal distribution is 
a statistically significant improvement over the single Gaussian 
at a confidence level of 89%. However, a much clearer result 
is obtained when histograms are built considering the bright- 



and the faint-SGBs separately (right panel of Fig. \T6[ : they 
have a different (averaged) C+N content, the faint-SGB hav- 
ing A(C+N)^7.64±0.24 and the bright-SGB A(C+N)^7.23±0.3 
(A(C+N)^7.80±0.19 and A(C+N)^7.47±0.26; using isochrone 
fitting temperatures, respectively, see Sect. 14. 3b . As an additional 
check, we performed a two-sample KS test computing the prob- 
ability that these two samples are drawn from the same parent 
distribution and found a pretty high {Pks ~0.03) significance of 
the difference in the bright-SGB and faint-SGB distribution of 
A(C+N). From Fig. [15] bright-SGB stars appear to be on aver- 
age more N-poor than then faint-SGB stars; assuming an N-O 
anticorrelation, the bright-SGB stars are more O-rich than the 
faint-SGB stars. Even though we do not provide oxygen abun- 
dances for our SGB stars, we can speculate on the C+N+O 
sum for the two SGB components assuming [O/Fe] values from 
available measurements. Various oxygen abundance determina- 
tions of RGB stars can be found in the literature. From Fig. [T3] 
we just note that large systematic differences exist between dif- 
ferent determinations of the O content (see also Tab. [3]i and 
we caution readers that assigning a reference [O/Fe] content to 
each SGB group could be naive at this stage. If we assume for 
the bright and faint component [O/Fe] ^0.1 dex and [O/Fe]^- 
0.2 dex, respectivel}[]3l the separation one sees in C+N con- 
tent virtually disappears when considering C+N+O. We found 
that the faint-SGB have A(C+N)^7.89+0. 14 and the bright-SGB 
A(C+N)^7.93±0.07. The distributio ns even swap when con- 
sidering O abundances suggested by lYong & Grundahll (l2008h 
([O/Fe]^0.5 dex and [O/Fe]^0.2 dex for the bright- and faint- 
SGB stars, respectivelj[3 We conclude that the bimodality we 
observe in the C+N sum does not necessarily imply or exclude 
a bimodality in the C+N+O content and more observations are 
needed to settle the case of NGC 1851. 

6. Summary and discussion 

We presented low-resolution spectroscopy for a large sample of 
MS and SGB stars in NGC 1851 with the goal of deriving C 
abundances (from the G band of CH) and N abundances (from 
the CN band at ~ 3883A) and investigating the chemical differ- 
ences between the two branches of the double SGB. We derived 
carbon and nitrogen abundances for 64 stars, whose spectra were 
obtained with FORS2 at VLT and IMACS at Magellan and ana- 
lyzed in a uniform manner. 

NGC 1851 is one of the most interesting GCs whose CMD 
displays a discrete structure at the level of the SGB and of 
the RGB. The photometric complexity is reflected in a peculiar 
chemical pattern that has only recently been investigated in de- 
tail. So far, the available abundance studies in NGC 1851 were 
limited to evolved stars that belong to the RGB (except for the 
study of iPancino et al This is the first time that a precise 

chemical tagging of C and N content is made for stars directly 
located in the bright- and faint-SGB component. The main re- 
sults of our analysis can be summarized as follows: 

- We derived CH and CN band index measurements for 
23 stars observed with IMACS, the spectrograph on the 

" w e de rived these averaged v alues from the works of lVillanova et al.l 
1 I2OIOI) and lCarretta et alj iioKt ) reported in Table[3| 

he re we no t e that our N abundances are systematically lower 
than Yong et al.l | |2009|) . in some case as much as 0.6-0.7 dex. We 
can attribute this discrepancy to (a) the different spectral resolution, 
(b) the different ev olutionary status of pro gram stars ( see F ig. 10 of 
iGratton et ani200(]|) and (c) the fact that in lYongetai] j2009l) N mea- 
surements come from the CN features at 8005 A. 
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Magellan I telescope . We added to our sample spectra from 
iPancino et aP (l2010l) . We were able to detect a large scat- 
ter and a hint of bimodality in the CN band strengths to- 
ward the brighter luminosities (refer to Fig. HJ. We did not 
report any clear anti-corellation from these index measure- 
ments (Fig.|5]). 

We performed spectral synthesis to separate the underlying 
C and N abundances from the CH and CN band strengths. 
Star-to-star strong variations with a significant range in A(C) 
and especially in A(N) were found at all luminosities from 
the MS (V ~ 20.1) up to the lower RGB (V ~ 18.6). C and N 
abundances are strongly anticorrelated, as would be expected 
from the presence of CN-cycle processing exposed material 
on the stellar surface (Fig.|7l). 

We us ed literature photometry in u, b, v, and y S tromgren 
bands dGrundahl et al.lll999l:ICalamida et al.ll2007h to define 
a new color index {{u-y)+{\'-b)). We found that the u versus 
{u-y) + {v-b) diagram is a powerful tool to identify the dou- 
ble RGB and SGB of NGC 1851 and showed that the faint- 
SGB is clearly connected with the red-RGB, while blue- 
RGB stars are l inked to the bright-SGB (Fig. [11] see also 
iHan et all l2009l) . Moreover, the relative frequency on the 
fainter/brighter SGBs (35% versus 65%) roughly matches 
the relative frequency of red-/blue-RGB stars selected in the 
u, {u-y) + (v - b) diagram (30% versus 70%; see Fig.fT2li. 
We photometrically defined blue- and red-RGB stars accord- 
ing to their position on this bimodal RGB sequence. We 
used i-elements, Na, O, and iron abundance that are avail- 
able from literature for some RGB stars of both populations 
to investigate their chemical content. The less populous red- 
RGB consists of Ba-rich La-rich stars and have, on aver- 
age, a higher Na abundance, while the bulk of Ba-poor La- 
poor stars belong to the blue-RGB. However, since we have 
demonstrated that the two RGB and SGB are photometri- 
cally connected, we can confidently extend these results to 
the two SGB components for these i-process elements not 
studied in this paper. 

Similarly, we isolated bona fide stars on the faint-SGB and 
bright-SGB using available photometry and analyzed their 
chemical composition. We noted a fully extended C-N anti- 
correlation for the bright-SGB stars, while faint-SGB stars 
tend to be richer in N, on average (Fig. [TsT l. The C-N pat- 
tern observed for SGB stars recalls the Na-O antic orrelation 
analyzed for RGB stars in previous papers. Specifically, the 
faint-SGB and bright-SGB samples are not completely su- 
perimposed on one another in the A(C), A(N) plane; but 
faint-SGB stars have, on average, a higher nitrogen abun- 
dance. This finding rules out the claims by ICarretta et alJ 
(1201 ih . who suggested that the faint-SGB is also nitrogen- 
poor. 

We analyzed the C+N sum for both bright-SGB and 
faint-SGB bona fide stars. Bright-SGB stars have 
A(C-i-N)^7.23+0.31 dex, while for the faint compo- 
nent A(C+N)^7.64.+0.24 dex. A diff^erence in log^(C+N) of 
0.4 dex as we find implies that the fainter SGB has about 2.5 
times the C-i-N cont ent of the brighter one. According to the 
ICassisi et a D ll2008h scenario, the faint-SGB is anticipated 
to have the higher CNO content. The current findings of 
increased C-hN content in the faint- SGB relative to the 
brighter one agree, in part, with the ICassisi eFal ] (I2008h 
results. However, we caution that the separation one sees in 
C+N content could significantly decrease or disappear when 
considering the Ch-N-i-O sum (as discussed in Sect. 15.3b . 



The general picture demonstrates that NGC 1851 shows an 
impressive resemblance to M 22. 

M 22 possesses a spread in s-process elements, iron content 
(although this is still debated for NGC 1851), and each of the 
two populations exhibits its own anticorrelation, with the ^-rich 
having on average higher C, N, and Na abundances. The chem- 
ical anomalies point to a bimodal SGB and RGB both for M22 
and NGC 1851. Similarly to NGC 1851, also for M22 the faint- 
SGB and the bright- SGB consist of i-rich and ^-poor stars (see 
iMarino et al ]l2012bh . 

Since the Na-O and the C-N anticorrelations alotie can be 
considered as the signature of multiple stellar populations, and 
both clusters are composed of two groups of stars with different 
i-element content (associated to the double SGB and RGB) pos- 
sibly with their own Na-O, C-N anticorrelations, we conclude 
that each group in turn is the product of multiple stellar forma- 
tion episodes. 

NGC 1851 and M 22 do not harbor only two stellar popula- 
tions (like normal GC) but have experienced a much more trou- 
bled star-formation history that resembles the case of wCentauri 
(see e. g. discussions in[Marino et a l. 2009; Da Costa et al. 200 ^ 
Da Costa & Marinoll201 iMRoederer et alJlMlllD' Antona et al.l 



201 Ih . 



genitor s ystems of u Cen, an d similarly M22 and NGC 1851 



ID' Antona et alJ (1201 ih suggested for a> Centauri a chemical 
evolutionary scenario where due to the large mass of the proto- 
cluster and its possible dark matter halo the material ejected by 
SNll may survive in a torus that co llapses back onto the clus- 
ter after the SN 11 epoch (see also iD'Ercole et alj|2008l) . The 
3D-hydro simulations by iMarcolini et alj ( l2006h show indeed 
that the collapse back includes the matter enriched by the SN 
11 ejecta. This scenario could be e asily extended to M22 and 
NGC 1851 (see lMarino et al.l2012ah . For u Centauri and M 22 it 
is tempting to speculate that enrichments in N and Na and deple- 
tion of C and O may have originated from the ejecta, collected 
in a cooling flow, of AGB stars that evolve in the cluster when 
the gas has been entirely exhausted by previous star-formation 

events. 

ID' Antona et alj ( 1201 Ih suggested that a poorly discussed site 
of s-nucleosynthesis that occurs in the carbon burning shells o f 
the tail of lower mass progenitors of SNll (e.g. lThe et al.ll200"7l) . 
may become particularly apparent in the evolution of the pro- 

ys , 

(see also lRoederer et ai1l201 ih . 

As an alternative possibility, NGC 1851 has been recently 
suggested t o be the merger-prod uct of two independent stellar 
aggregates (Ivan den Berghl[l996l) . While this possibility seems 
unlikely for globular clusters in the Galactic halo, an origin as a 
merger product of two independent star clusters cannot be ex- 
cluded in dwarf gala xies. In this case, numerical simulations 
(iBekki & Yong|l20rih showed that two clusters can merge and 
form the nuclear star cluster of a dwarf galaxy. After the parent 
dwarf galaxy is accreted by the Milky Way, its dark matter halo 
and stellar envelope can be stripped by the Galactic tidal field, 
leaving behind the n ucleus (i.e., NGC 1851 ) and a diffuse stellar 

halo (as observed bv lOlszewski et alll2009l) . 

As already mentioned in the introduction, ICarretta et alJ 
(l20Tll) associated the i-rich and the i-poor populations to the 
bright-SGB and the faint-SGB, respectively, with the bright- 
SGB having also higher N and Na abundances. According to 
Carretta and collaborators, the possibility that the faint-SGB is 
CNO enhanced should be excluded, demonstrating that the split 
is caused by an age difference of ~ 1 Gyr between the two pop- 
ulations. In this paper we have shown instead that the faint-SGB 
is made of N-rich and probably ^-rich stars and bright-SGB 
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stars are N-poor and probably 5-poor. While we added impor- 
tant pieces of information to the general picture, our results do 
not provide a conclusive answer on the occurrence of a merger 
in NGC 1851 and suggest that the measurement of the overall 
C+N+O abundance as well as a precise determination of the spa- 
tial distribution of the multiple SGBs and RGBs are still manda- 
tory to shed light on the star-formation history of this GC. 
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Lardo et al.: C and N abundances of MS and SGB stars in NGC 1851. 
Table 1. Index measurements for NGC 1851 stars 



ID Ra Dec V (V-I) CN errcN 5 CN CH errcH (5CH 

(deg) (deg) (mag) (mag) (mag) (mag) (mag) 



1 1219 


78 


.5334822 


-40. 


AO n A A cc\ 

,02 /4459 


18, 


.643 


0. 


837 


1 1755 


78. 


C Af\A AO 


-40. 


A1 00/;0/C 

,0328686 


18, 


.774 


0. 


818 


12485 


78. 


.5406833 


-40. 


,0253661 


19 


.001 


0. 


697 


12925 


78 


,344 /oD / 


-40. 


AO T /I O C /I 

,0234854 


19 


.297 


0. 


660 


13062 


78 


C C 1 TAT7 


-40. 


AQ A 1 1 AC 

.0301395 


18, 


.624 


0. 


827 


13872 


78. 


C COT /I 


-40. 


f\A A C "2 AT 

,044530 / 


19 


.263 


0. 


655 


15182 


78. 


C/^A"? Ano 

.5603478 


-40. 


f\A £^C\C O O 

.0469688 


18, 


.646 


0. 


857 


15490 


78 


.5599836 


-40. 


f\A OOCC\ A 

,0488694 


18, 


.977 


0. 


,714 


16047 


78 


,5459865 


-40. 


AC O 1 AAO 

,0581902 


18, 


.610 


0. 


844 


21)293 


16. 




A n 
-41). 


,0632061 


1 o 

Is, 


.OiU 


A 

0. 


OQ O 

03S 


40017 


78. 


.4595330 


-40. 


1 CO COAT 

,152529 / 


19 


.695 


0. 


615 


4U02U 


la 


/I TO /I /icn 
.4/24050 


A n 


1 C 1 1 A 1 A 

,151 1010 


1 o 

lo. 


O /I 1 

.041 


A 

0. 


Tn^ 
/96 


40022 


78 


.4332081 


-40. 


1 C A /I T/; A 

,1504760 


19 


.888 


0. 


,652 


A nmo 
4U025 


TO 

/O. 


.4o954Ud 


A n 
-41). 


1 AC\A ^^CO 

, 1494668 


1 A 

19 


ACiA 

.494 


A 

0. 


63J 


40051 


78. 


.4631667 


-40. 


1 A 'JC\ A A A 

,1439444 


18, 


.792 


0. 


,796 


4U062 


no 

la. 


AQi A A on 
.48J448 / 


A n 


A A A A r\A n 
,141404/ 


1 A 

19 


1 ^ 1 

. lot 


A 

0. 


664 


40072 


78 


An Of\C A1 

,4780543 


-40. 


1 Q AT C/; O 

1393568 


19 


.934 


0. 


648 


A nmo 
41)0/ 6 


TO 

16. 


,4jDo / lo 


A n 
-41). 


1 Q o 1 i^cn 
130166/ 


1 A 

19 


Ai^A 

.4o4 


A 

0. 


029 


40083 


78. 


,430/625 


-40. 


1 Q T 1 C O /I 

13 / 1 584 


19 


.305 


0. 


647 


A nnoo 
4l)005 


TO 

16. 


,4oy!)o40 


A n 
-41). 


1 Q cm/1 
1359240 


1 A 

19 


n AO 
.240 


A 

0. 


f^An 
,04/ 


A r\r\r\ a 
40094 


78 


A nn A 1 on 
,4/ /43o / 


-40. 


1 Q /I TT OT 

134//8/ 


19 


.490 


0. 


,640 


40097 


78 


.4733785 


-40. 


1341169 


20, 


.053 


0. 


658 


40100 


78. 


,44v5d1 1 


-40. 


1 Q /I A AAA 

1 340000 


18, 


.924 


0. 


,751 


40117 


78. 


A COAT £^ C 

.4580365 


-40. 


1314412 


19 


.358 


0. 


652 


40123 


78 


A £iC /I TCA 

.4664750 


-40. 


13061 1 1 


19 


.006 


0. 


,709 


40133 


78 


A A A A ACiC 

,4144050 


-40. 


1 O OTC AA 

128 /500 


19 


.193 


0. 


667 


40153 


78. 


A ^1 A mo 
,45349/2 


-40. 


1264167 


19 


.479 


0. 


632 


40167 


78. 


A A Cn d^C\A 

,445 /d94 


-40. 


1245556 


18, 


.932 


0. 


,706 


40186 


78 


,404 /8oi 


-40. 


1218056 


18, 


.869 


0. 


,817 


40191 


78 


,4/022o3 


-40. 


1 O 1 OAA/; 

,1212996 


19 


.598 


0. 


,617 


40196 


78. 


A ATjmn 

.4923929 


-40. 


1 OA /I AAC 

1204995 


19, 


.649 


0. 


622 


40197 


78. 


/I T /I C /I OA 

.4345480 


-40. 


1 OA/1 AOC\ 

,1204489 


19 


.412 


0. 


657 


40235 


78 


C 1 Al AC\ 

.5103649 


-40. 


1 167899 


19 


.921 


0. 


,646 


40239 


78 


A c\ocno A 

.4986784 


-40. 


1 165524 


19 


.562 


0. 


,620 


40241 


78. 


A A 'J AO C A 

,4430254 


-40. 


11652/2 


19 


.693 


0. 


622 


40247 


78. 


A OOO CA 

.4682250 


-40. 


1159/22 


18, 


.705 


0. 


837 


40271 


78 


.5131639 


-40. 


11/1 A AAA 

1 140000 


18, 


.866 


0. 


,844 


40303 


78 


,44i 130d 


-40. 


1118056 


19 


.033 


0. 


686 


40340 


78. 


/m A 1 r\A 
,4/30104 


-40. 


1 AAT A 1 A 

1093919 


19 


.913 


0. 


642 


40344 


78. 


C AO OAA 1 

.5028991 


-40. 


1 AAAC C A 

,1090554 


19 


.438 


0. 


633 


40348 


78 


,445680/ 


-40. 


1 AOA A on 
,108948/ 


19 


.820 


0. 


646 


40376 


78 


A nnr\nf\A 
,4/ /0/04 


-40. 


, 106/ /91 


19 


.930 


0. 


659 


40378 


78. 


CO C AC O T 

.5259583 


-40. 


1 ATAAAA 

,10/0000 


18, 


.740 


0. 


825 


40385 


78. 


C A A 1 C Q O 

,5091532 


-40. 


1 A/; 1 /;OA 

,1061629 


19 


.806 


0. 


630 


Af\A^ A 

40424 


78 


.4685664 


-40. 


1 Al /I C O A 

1034580 


19 


.438 


0. 


,657 


40431 


78 


.5162111 


-40. 


1 Al 1 1 OA 

,1031389 


18, 


.881 


0. 


811 


40465 


78. 


An A iL^n o 
.4/462/8 


-40. 


1 A 1 T T TI 

1013333 


18, 


.976 


0. 


,757 


40504 


78. 


/I O C /I C 1 

.4685451 


-40. 


AAAO A AA 

,0992400 


19 


.454 


0. 


597 


40507 


78 


A TA/^TOO 

,4/06/22 


-40. 


AAAQ t 1 

099361 1 


19 


.173 


0. 


607 


40508 


78 


C AA A Q Al 

.5094543 


-40. 


AAAAO/^ 1 

0990863 


19 


.602 


0. 


625 


40545 


78. 


A n cnnn o 

A/di / /8 


-40. 


AAT COTO 

,09 /52/8 


18, 


.985 


0. 


,728 


40571 


78. 


c A /m CO 
,504/352 


-40. 


aa/;aa a o 
,0960942 


19 


.437 


0. 


611 


40575 


78, 


A /'AO A f\'~l 

.4698407 


-40. 


AA/'A 1 A/' 

0960196 


20, 


.013 


0. 


,649 


40620 


78, 


.4532694 


-40. 


.0945000 


18, 


.640 


0. 


,835 


40665 


78, 


.5128445 


-40. 


.0926724 


18, 


.995 


0. 


,712 


40679 


78, 


.4762083 


-40. 


.0924444 


18, 


.716 


0. 


,803 


40709 


78, 


.4309537 


-40. 


.0912317 


19 


.166 


0. 


,657 


40715 


78, 


.4849583 


-40. 


0910365 


19, 


.908 


0. 


,644 


40756 


78, 


.4788241 


-40. 


.0895444 


19 


.280 


0. 


,647 


40827 


78, 


.5001750 


-40. 


.0876944 


18, 


.762 


0. 


,836 


40863 


78, 


.5070126 


-40. 


.0863515 


19 


.393 


0. 


608 


40874 


78, 


.4722925 


-40. 


.0860825 


19 


.455 


0. 


616 


40919 


78, 


.5234044 


-40. 


.0849874 


19 


.496 


0. 


,637 


40978 


78, 


.4762989 


-40. 


.0829703 


19 


.919 


0. 


,640 



A Q O 1 

-U.381 


A A/1 T 
U.U4 / 


A AA^ 

-0.006 


A 


O AO 

.8U8 


A A/1 ^ 

U.U40 


A AT O 
0.030 


-0.237 


0.058 


0.134 





.738 


0.054 


-0.038 


A 1 

-0.2ol 


A A^Q 

U.Uo3 


A 1 A^ 

U. lOo 


A 



TT /I 

. /34 


A AiiO 

U.U62 


A A/1 O 

-0.042 


-0.380 


0.056 


-0.002 


0, 


.781 


0.039 


0.017 


A QT 1 
-U.J / i 


U.0o3 


A A A /I 

-0.004 


A 


O 1 ^ 

.816 


A ACC 

U.U53 


A ATA 

0.039 


A Q OO 

-U. j8o 


A AQ O 

U.U3o 


A A 1 ^ 

-U.Uio 


A 

U, 


OTO 

.8 16 


A A/1 T 

U.U43 


A 1 AO 

0.102 


A Q 1 /I 


A A/1 

U.U40 


A t\CA 

0.064 


A 



nin 


A AOA 

U.U29 


A AO£ 

-0.020 


A TAT 

-0.2U/ 


A AT 1 
U.U31 


A 1 T 1 
U. 1 / 1 


A 



O 1 A 
■ OlU 


A A/1A 


A A/1 ^ 

0.040 


A /1A/C 

-0.4U6 


A AT 1 

0.023 


A AOA 
0.021) 


A 

0, 


.6/9 


A A/1 T 

U.U43 


A AT T 

-0.033 


-0.230 


A AO A 

U.U24 


A 1 A^ 

U. lOo 


A 

0, 


O /I 1 

.041 


A A/1 C 

U.U43 


A AT^ 

0.030 


-0.336 


0.023 


0.105 





.669 


0.042 


-0.053 


-0.418 


0.023 


-0.007 


0, 


.716 


0.036 


0.002 


A Q CT 

-0.35 / 


A ACT 

0.052 


A A 1 O 

O.Olo 


A 



oil 
■ Oil 


A A/1 T 

U.U43 


A A/1 A 

0.041) 


-0.324 


0.023 


0.062 


0, 


.772 


0.036 


0.028 


-0.412 


0.025 


0.032 


0, 


.684 


0.049 


-0.042 


A /I 1 Q 


A AC 1 

U.U3 1 


A A/1 A 

-0.049 


A 


TOC 

. /23 


A AC^ 

U.U56 


A A/1 O 

-0.040 


A QAA 

-U.390 


A AT 1 

U.U33 


A AAT 
0.00/ 


A 

0, 


cn A 
.0/4 


A A/1 C 

U.U43 


A ACT 

-0.053 


A QAC 

-U.395 


A AO A 

U.U24 


A AAT 

-0.003 


A 

U, 


/COA 

.6oU 


A A/1 O 

U.U42 


A ACT 

-0.053 


-0.464 


0.032 


-0.053 


0, 


.693 


0.045 


-0.021 


-0.434 


0.025 


0.019 


0, 


.755 


0.050 


0.015 


A QAA 

-0.390 


A A£ 1 

O.Ool 


A A 1 T 

-0.01 / 


A 

U, 


TAA 
. /99 


A ATA 
U.U39 


A AO A 

0.024 


-0.345 


0.026 


0.056 





.720 


0.049 


-0.002 


-0.387 


0.054 


-0.016 


0, 


.716 


0.042 


-0.060 


-0.380 


0.036 


0.039 


0, 


.699 


0.048 


-0.012 


-0.436 


0.022 


-0.013 


0, 


.662 


0.040 


-0.049 


-0.355 


0.020 


0.050 





.733 


0.031 


0.015 


-0.272 


0.021 


0.171 


0, 


.723 


0.031 


-0.002 


-0.415 


0.019 


0.001 


0, 


.651 


0.036 


-0.061 


A /ino 
-0.4U0 


A AO 1 

U.U23 


A A 1 O 

O.Olo 


A 




cnz 
.6/5 


A A/1A 
U.U4U 


A ATT 

-0.03 / 


-0.375 


0.039 


0.001 





.814 


0.058 


0.046 


A T /I O 

-U.240 


A ACC 


A 1 O^ 

U. 12o 


A 

0, 


AOC 

.923 


A A/1A 
U.U4U 


A 1 CO 

0.132 


A /1 1 n 
-0.419 


A A 1 n 
0.019 


A AOT 

0.023 


A 

0, 


TO/1 

. /24 


A ATT 
U.IJ3 / 


A AAA 

-0.000 


-0.444 


0.027 


-0.037 


0, 


.726 


0.037 


0.010 


A QT 1 

-U.3 / 1 


A AOT 
U.U2 / 


A t\CA 

0.064 


A 




TO/1 

. /24 


A ATT 
U.U3 / 


A AAT 
0.00/ 


-0.401 


0.034 


0.043 


0, 


.754 


0.055 


0.028 


-0.449 


0.020 


-0.015 


0, 


.735 


0.044 


0.019 


A QAn 
-U.3U9 


A AOT 
U.U2 / 


A AAO 

0.090 


A 




T/1 A 
. /4U 


A ATT 
U.U3 / 


A AO A 

0.024 


-0.476 


0.101 


-0.104 


0, 


.733 


0.156 


-0.043 


-0.405 


0.020 


0.003 





.724 


0.044 


0.008 


A QA/i 

-U.396 


A AQ/i 

U.U36 


A AOT 

-0.02/ 


A 

u. 


TCO 

. /32 


A AC 1 


A AOC 

-0.023 


-0.350 


0.025 


0.069 





.689 


0.038 


-0.022 


-0.380 


0.035 


-0.008 


0, 


.765 


0.043 


-0.011 


-0.321 


0.027 


0.086 





.707 


0.039 


-0.010 


-0.498 


0.025 


-0.048 


0, 


.740 


0.042 


0.005 


-0.458 


0.030 


-0.081 


0, 


.870 


0.042 


0.105 


-0.443 


0.029 


-0.070 


0, 


.765 


0.037 


-0.007 


-0.250 


0.038 


0.126 





.821 


0.039 


0.053 


-0.494 


0.020 


-0.108 


0, 


.707 


0.044 


-0.036 


-0.388 


0.025 


0.054 





.700 


0.038 


-0.024 


-0.427 


0.036 


-0.032 


0, 


.730 


0.040 


0.001 


-0.396 


0.018 


0.007 





.737 


0.040 


0.018 


-0.341 


0.027 


0.067 





.721 


0.039 


0.005 


-0.424 


0.023 


0.019 





.691 


0.047 


-0.034 
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Ka 


Dec 


V 


(V-ij 




errcN 


L-JN 


Lri 


errcH 


Crl 






(deg) 


(deg) 










(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


41003 


78 


AC\A£^OA A 

,494oo44 


-40.082507 1 


1 A 

19, 


ACA 

.050 


0.729 


-0.495 


0.029 


-0.118 




0.037 


0.005 


41018 






A(\ nooQO<;;n 
-41). 0823209 


1 

18, 


occ 
.853 


A 1 A 

0.819 


A Q/;c 
-0.3d5 


A ATT 
0.03/ 


A AAA 
0.009 


A 1 
O.8IO 


A A/1 1 

0.041 


A A/1 1 

0.043 


41 108 


78. 


A C O "2 CO 

,4583652 


Af\ AOAATCO 

-40.0800/52 


1 A 

19 


nzA 
. /D4 


0.656 


-0.423 


0.025 


0.007 


0.722 


0.038 


0.008 


41 loD 


TO 

la. 


,jUU945i 


A(\ ATOCOTO 

-40.0/858 /8 


1 A 

19 


1 1 A 
. 1 lU 


r\ CAC 
O.D4o 


A /I 1 

-U.418 


A AT A 
U.U24 


A AQ^ 

-0.030 


A TAA 

0. Iw 


A AQ C 

0.033 


A A/1 Q 

-0.043 


41213 


78 


c 1 /; yi n /lo 

.5104942 


Af\ AT0T110 

-40.0/82312 


1 A 

19 


1 ZA 

.134 


0.695 


-0.371 


0.024 


-0.002 


0.725 


0.038 


-0.052 


412/y 


TO 

/5. 


,42dUd92 


A(\ AT^OOQn 

-40.0/OZ2J9 


1 A 

19 


TTT 

.111 


O.DjJ 


A Q A 1 

-U.391 


A AT/; 
0.U2O 


A AAQ 

0.003 


A /iOQ 

0.083 


A AQA 

0.039 


A A/IT 

-0.04/ 


41325 


78. 


.4670381 


/1A ATCO/;C/1 

-40.0/52o54 


1 A 

19 


A OA 

.48U 


0.636 


-0.360 


0.027 


0.050 


0.702 


0.040 


-0.012 


41 JjU 


TO 

la. 


.4/ j9yD4 


/I A AT/I /iTC /I 

-40.0/4D/D4 


1 A 

19 


Q AO 

.308 


0.D43 


A CAA 


A A 1 

U.U18 


A 1 1 T 

-0. 1 12 


0.003 


A AQA 

0.039 


A AiiQ 

-0.0o3 


41372 


78 


.4922083 


A A AT /I Q t 1 

-40.0/4361 1 


1 

18, 


.III 


0.856 














A 1 ceo 


TO 

la. 


/I TnooTn 
.4/028/y 


A A AT 1 AT 1 T 

-4U.U/10/1 / 


1 A 

19 


TO 1 

. /21 


O.OJO 


A Q 00 

-U.388 


A ATQ 

0.023 


A A/1 A 

0.040 


A <;;oA 
0.D8U 


A A/IT 

0.04/ 


A AQ 1 

-0.033 


41610 


78. 


.454O060 


A A AT AA/; 1 A 

-40.0/00619 


1 A 

19 


1 in 

.HI 


0.684 


-0.454 


0.036 


-0.070 


0.743 


0.040 


-0.004 


A 1 f^ClA 

41d94 


TO 

la. 


.5000889 


/1A A/iOTOTO 

-40.008 1 111 


18, 


.736 


A /I T 

0.84/ 














4.1 807 


78 


4.07 1 (S7Q 




18, 


.792 


U.OJO 


-n ^67 


0.025 




yt.y 






41835 


78. 


.4320835 


-40.0662831 


19 


.012 


0.707 


-0.345 


0.018 


0.029 


0.700 


0.040 


-0.068 


41884 


78. 


.5039348 


-40.0656429 


18, 


.978 


0.775 














42073 


78, 


.4802671 


-40.0628473 


19 


.304 


0.640 


-0.418 


0.021 


-0.022 


0.737 


0.038 


0.010 


42195 


78, 


.4415574 


-40.0610284 


19 


.622 


0.643 


-0.394 


0.024 


0.026 


0.723 


0.035 


0.012 


42623 


78, 


.4582694 


-40.0551667 


18, 


.971 


0.734 


-0.413 


0.031 


-0.041 


0.715 


0.049 


-0.060 


42785 


78, 


.4376004 


-40.0528185 


19 


.496 


0.616 


-0.418 


0.035 


-0.007 


0.722 


0.043 


0.008 


42865 


78, 


.4964339 


-40.0519295 


18, 


.869 


0.794 














43014 


78, 


.4127639 


-40.0498611 


19 


.197 


0.623 


-0.379 


0.047 


-0.011 


0.722 


0.056 


-0.055 
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Table 2. Atmospheric parameters and carbon and nitrogen abundances for NGC 1851 stars 



ID 


Teff 

(K) 




A(C) 


eA(C) 


A(N) 


eA(N) 


11219 


5193 + 90 


3.4 


6.42 


0.14 


7.29 


0.31 


11755 


5247 + 93 


3.5 


6.51 


0.12 


7.15 


0.34 


12485 


5620 ± 116 


3.8 


6.64 


0.13 


7.58 


0.53 


12925 


5746 ± 124 


3.9 


6.88 


0.14 


7.73 


0.34 


13062 


5221 ± 92 


3.4 


6.14 


0.14 


7.50 


0.30 


13872 


5764 ± 125 


3.9 


7.10 


0.18 


7.17 


0.35 


15182 


5139 ± 87 


3.4 


6.20 


0.14 


7.45 


0.23 


15490 


5563 ± 112 


3.7 


7.01 


0.13 


6.36 


0.30 


16047 


5174 + 90 


3.4 


5.90 


0.14 


7.69 


0.26 


20295 


5191 + 90 


3.4 


6.25 


0.12 


7.48 


0.26 


40017 


5909 ± 135 


4.1 


6.55 


0.13 


7.12 


0.26 


40020 


5310 + 97 


3.6 


6.47 


0.12 


7.12 


0.19 


40022 


5774 ± 126 


4.2 


5.47 


0.14 


8.67 


0.26 


40062 


5732 ± 123 


3.9 


6.57 


0.11 


7.82 


0.24 


40072 


5789 ± 128 


4.2 


6.44 


0.14 


7.44 


0.26 


40083 


5792 ± 127 


3.9 


6.46 


0.12 


7.20 


0.25 


40088 


5792 ± 127 


3.9 


6.40 


0.12 


7.52 


0.25 


40094 


5817 ± 129 


4.0 


6.32 


0.12 


7.44 


0.23 


40097 


5753 ± 125 


4.2 


6.61 


0.14 


7.34 


0.26 


40117 


5774 ± 126 


4.0 


6.66 


0.12 


7.37 


0.25 


40123 


5580 ± 102 


3.7 


6.60 


0.15 


7.43 


0.53 


40191 


5902 ± 135 


4.1 


6.09 


0.13 


8.00 


0.27 


40197 


5757 ± 125 


4.0 


6.33 


0.13 


7.99 


0.24 


40235 


5796 ± 128 


4.2 


6.42 


0.14 


8.07 


0.24 


40239 


5890 ± 135 


4.1 


6.22 


0.13 


7.66 


0.25 


40340 


5810 ± 129 


4.2 


6.59 


0.13 


7.89 


0.22 


40344 


5843 ± 130 


4.0 


6.64 


0.12 


7.28 


0.27 


40348 


5796 ± 128 


4.2 


6.54 


0.13 


6.32 


0.24 


40376 


5750 ± 124 


4.2 


6.43 


0.15 


7.52 


0.25 


40385 


5854 ± 132 


4.2 


6.56 


0.12 


7.01 


0.26 


40424 


5757 ± 125 


4.0 


6.45 


0.13 


7.40 


0.24 


40504 


5977 ± 140 


4.1 


6.76 


0.13 


7.42 


0.26 


40507 


5939 ± 117 


3.9 


6.72 


0.16 


7.85 


0.29 


40508 


5872 ± 133 


4.1 


6.25 


0.14 


7.47 


0.30 


40545 


5518 ± 100 


3.7 


6.41 


0.16 


7.64 


0.28 


40571 


5924 ± 137 


4.0 


6.33 


0.13 


7.99 


0.25 


40575 


5785 ± 127 


4.2 


6.61 


0.15 


6.66 


0.30 


40620 


5199 + 86 


3.4 


6.57 


0.14 


6.71 


0.23 


40665 


5570 ± 113 


3.7 


6.50 


0.12 


7.05 


0.22 


40709 


5757 ± 125 


3.9 


6.62 


0.13 


7.95 


0.25 


40715 


5803 ± 129 


4.2 


6.51 


0.13 


7.41 


0.26 


40756 


5792 ± 127 


3.9 


6.48 


0.12 


7.71 


0.24 


40863 


5935 + 130 


4.0 


6.31 


0.17 


7.58 


0.36 


40874 


5905 ± 135 


4.0 


6.69 


0.14 


6.96 


0.27 


40978 


5817 ± 130 


4.2 


6.50 


0.18 


7.41 


0.29 


41003 


5515 ± 109 


3.7 


6.47 


0.13 


6.81 


0.22 


41018 


5244 + 93 


3.6 


6.41 


0.12 


7.05 


0.26 


41108 


5760 ± 125 


4.1 


6.51 


0.14 


7.47 


0.26 


41185 


5796 ± 127 


3.9 


5.78 


0.13 


7.52 


0.27 


41213 


5626 ± 116 


3.8 


6.44 


0.13 


7.77 


0.25 


41279 


5843 ± 131 


4.0 


6.36 


0.14 


7.60 


0.29 


41325 


5832 ± 130 


4.0 


6.66 


0.14 


7.15 


0.25 


41350 


5806 ± 128 


4.0 


5.49 


0.14 


8.57 


0.25 


41558 


5832 ± 130 


4.1 


6.54 


0.14 


7.01 


0.28 


41610 


5663 ± 119 


3.8 


6.45 


0.13 


6.81 


0.22 


41807 


5136 + 87 


3.5 


6.67 


0.12 


6.37 


0.21 


41835 


5586 ± 114 


3.8 


6.45 


0.13 


6.63 


0.23 


41884 


5372 + 97 


3.7 


6.34 


0.31 


7.42 


0.36 


42073 


5817 ± 129 


4.0 


6.37 


0.13 


7.72 


0.25 


42195 


5806 ± 128 


4.1 


6.63 


0.13 


7.11 


0.27 


42623 


5499 ± 99 


3.7 


6.94 


0.13 


6.98 


0.22 


42785 


5905 ± 130 


4.1 


6.50 


0.14 


7.80 


0.26 


42865 


5316 ±93 


3.6 


6.56 


0.20 


7.00 


0.30 


43014 


5879 ± 115 


3.9 


6.86 


0.18 


7.78 


0.54 
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